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Equilibrium canonical distribution in statistical mechanics assumes weak system-bath coupling 
(SBC). In real physical situations this assumption can be invalid and equilibrium quantum statistics 
of the system may be non-canonical. By exploiting both polaron transformation and perturbation 
theory in a spin-boson model, an analytical treatment is advocated to study non-canonical statistics 
of a two-level system at arbitrary temperature and for arbitrary SBC strength, yielding theoretical 
results in agreement with exact Monte-Carlo simulations. In particular, the eigen-representation of 
system's reduced density matrix is used to quantify non-canonical statistics as well as the quantum- 
ness of the open system. For example, it is found that irrespective of SBC strength, non-canonical 
statistics enhances as temperature decreases but vanishes at high temperature. 

PACS numbers: 03.65.Yz; 03.67.Mn; 05.30.-d 



I. INTRODUCTION 

Equilibrium canonical distribution is a fundamental 
result in statistical mechanics, but with an implicit as- 
sumption that the system-bath coupling (SBC) is van- 
ishingly weak. In real physical situations, such as light- 
harvesting systems [l [, super-conducting qubits [2[ , and 
atom-cavity systems 3|, this weak SBC assumption can 
be invaUd. Then a good separation between the system 
and the bath is lost, leading to non-canonical equilib- 
rium statistics for the system (though the system plus 
the bath as a whole still satisfies a canonical distribu- 
tion). At present non-canonical statistics is a subject of 
great interest. It challenges intuitive concepts in statisti- 
cal mechanics Q and brings about new understandings of 
open quantum systems. Because canonical statistics can 
be recovered in the classical limit for a wide class of mi- 
croscopic open-system models non-canonical statis- 
tics can also serve as an indicator of the quantumness of 
open systems. 

To better understand non-canonical statistics and for 
future bench-marking purposes, an analytical treatment 
of non-canonical statistics at arbitrary temperature and 
also for the whole range of SBC strength is important. 
This task is achieved here in a spin-boson model by ex- 
ploiting both the polaron transformation technique and 
perturbation theory, with results applicable to almost the 
entire parameter space (certain conditions to be elabo- 
rated below). Though our theory does provide all the 
information of system's reduced density matrix (RDM), 
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one simple measure is employed to quantify the degree 
of non-canonical statistics. Specifically, we use the angle 
by which the diagonal representation of RDM is rotated 
away from the energy eigen-representation of the system 
as the measure. Evidently, this angle is zero for a canon- 
ical distribution because the associated RDM is diagonal 
in energy basis states. It is then shown how an increas- 
ing SBC strength or a decreasing temperature strength- 
ens non-canonical statistics. Interestingly, at sufficiently 
high temperature, our theory reveals that the standard 
canonical statistics will be recovered asymptotically. The 
theoretical results are found to be in full agreement with 
exact Monte-Carlo simulations. The observation that 
non-canonical statistics vanishes in the high temperature 
regime resonates with the view that non-canonical statis- 
tics reflects the quantumness of the open system, which is 
expected to diminish at high temperature. Furthermore, 
we exploit our simple non-canonical statistics measure to 
identify an interesting temperature scale, at which the 
non-canonical equilibrium statistics is most sensitive to 
temperature variations. 



II. THEORY 

The total Hamiltonian Ht of a system plus a bath is 
written as 
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where the three terms represent Hamiltonians of the sys- 
tem, the bath, and SBC, respectively. The system and 
the bath as a whole is assumed to be in contact with 
a super-reservoir at temperature T. Long after quan- 
tum relaxation takes place, the system and the bath as 
a whole reaches a canonical equilibrium state pt- The 
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equilibrium RDM of the system is given by 

Ps=tTB[pt]= ^^^^-,H,y (2) 

where /3 = ^'"[■] denotes the trace over the system 

as well as the bath, and trB[-] denotes the trace over 
the bath only. Due to the non-vanishing contribution of 
HsB to Ht, in general the equilibrium RDM is no longer 
a canonical distribution of the system alone It is 
hence of importance to develop a theoretical treatment 
to systematically investigate the dependence of ps upon 
temperature T and the magnitude of Hsb- We achieve 
this task here for a spin-boson model with the Hamilto- 
nian (we set ?i = 1) 
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where ai {i = x, y, z) are the usual Pauli matrices, e is 
the energy splitting between two levels of the system, 
and A is the tunneling matrix element. The bath is 
modeled as a set of harmonic oscillators with frequencies 
LOk, and the couplings to the spin system are denoted by 
gk- The properties of the bath are fully characterized by 
its spectral density, namely, J{uj) = T^J^k 9k^(^ ^ '^k)- 
For a fully analytical development we first assume a 
super-ohmic spectral density with exponential cut-off, 
i.e., J(w) = 7CJ'^e^'^/"'=, where 7 characterizes the SBC 
strength and Uc is the cut-off frequency of the bath. 
Though such an open-system model is standard and (de- 
ceptively) simple, even its stationary properties do not 
have exact solutions (not to mention its dynamics). In- 
deed, only the case of a single-mode spin-boson model 
with e = can be regarded as integrable |8| . As such, per- 
turbation theory becomes one of the few options. Note, 
however, a naive perturbation theory would fail as our 
central concern is beyond the weak SBC regime. 

To capture the effects of a finite Hsb, we exploit a 
standard paloron transformation as our starting step. 
To that end, we transform the total Hamiltonian Ht to 
Ht = e^Hte~^ (tildes denote operators in the trans- 
formed polaron picture) where 



F = a, 



i^k 



^-(bi-bk). 
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We then obtain (up to a constant) 
Ht = Hs + Hb + Hsb , 



with Hs = fcTz + -^f^x, Hb = Hb, and Hsb = cr^K + 
<JyVy. Two remarks are in order: (i) for Hs, the <Tx 
(tunneling) term is now renormalized by SBC, because 



Ai^/A = i? = exp 
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du J(lo) 
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Using the super-ohmic J(w) mentioned above, we obtain 
R = cxp[-^(2i/''(l/^Wc) - t^lP^)], with V' being the 
derivative of the digamma function; and (ii) Hsb, the 
SBC in the paloron picture, assumes a form very different 
from Hsb, because the bath operators entering into SBC 
are now 



Vx 



{D\ + Dl - 2R), 
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(7) 
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withD±^exp[±^{bl~bk)] and {DI)hs = i? @. We 
stress that so far our procedure is exact. As seen below, 
this polaron transformation is advantageous in treating 
cases of non-vanishing SBC strength because the corre- 
lation functions of Vx and of Vy decay faster for a larger 
7. For later use we also introduce density matrices in 
the polaron picture, i.e., p = e^pe~^, where p can be 
the RDM Ps or the full density matrix pt- For example, 



Pt 



^ and ps = tTB[pt]- 



tr[e-'9«t] ■ 

We now proceed with the calculations of the matrix el- 
ements of Ps, expressed in terms of the Uz basis. The di- 
agonal elements p^g ^'^d arc given by ^^{l^iTs^^zPs]) 
and the off-diagonal element p^g is given by lisV'-Ps] ~ 
tr [cT-pt]. Noticing that ii s[<yzPs] = trsio-^Ps], one may 
obtain the diagonal elements directly if ps, the RDM in 
the polaron picture, is analytically known. On the other 
hand, because (t_ operator does not commute with the 
polaron transformation operator F defined above, calcu- 
lating p^g in the polaron picture is more involving. Nev- 
ertheless, we find 



pf = tr[cr_pt] = tr[o-_pt] = tr[o-- 



-Dipt] 



(9) 



indicating that p^g can be still obtained from pt , but not 
from Ps- 

Having observed how all RDM elements may be eval- 
uated in the polaron picture, we turn to the central ex- 
pression e~^^' in Pt. We treat /3 as an imaginary time 
and exploit a "time" -dependent perturbation theory in 
terms of Hsb- We first briefly mention the calculations 
of the RDM diagonal elements, a subject of a recent tech- 
nical study by us p^ . In particular, due to the fact that 
irB[HsB\ ~ 0, the leading-order contribution of Hsb is of 
the second order, from which we have ps ~ Ps,(o)~^ Ps,{2)^ 
where [ill 
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trs[c-^^^] and Zs,(2) = tis[A]. In addition, by 



the operators in imaginary time are defined as 0(t) 
e^^oOe-^^n, where i?o = Hs+Hb- The bath correlation 
functions C„(t) = {Vn{T)Vn)HB are given by 



2), 
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with </>(r) =4/„°' 



A straight- 



poo du> cosh[|(ff-2T)c 
1^2 sinh(,9w/2) 

forward inspection of these expressions (in particular, the 
term A^e'^^^^) shows that as 7 (the SBC strength) in- 
creases, the correlation functions Cn(r) always decreases 
exponentially with 7, thus enhancing the perturbation 
theory in imaginary time. Indeed, in the strong SBC cou- 
pling limit, the second-order correction pg approaches 
zero. This roughly explains how our polaron-picture- 
based perturbation theory, by construction, may be suit- 
able for treating finite SBC strength. It is now also clear 
that our theory is pcrturbative in terms of Hs b , but non- 
perturbative in terms of Hsb- We also note that though 
the integral expression in Eq. ()lip is complicated, it in- 
volves 2x2 matrices only. The final expression for the 
trace trs\aj.ps] is hence also analytical. Detailed calcu- 
lations [lO| show that the RDM diagonal elements thus 
obtained are accurate so long as the tunneling element 
A is not large as compared to the bath cut-off frequency 
ujc (therefore not a slow bath). 

The crucial task in theory here is to explicitly evaluate 
the off-diagonal element of ps via tr[o'__D^/5t]. In this 
case, due to the correlation between the system and the 
bath, the first-order contribution of Hsb to e^^^* and 
hence to pt are already nonzero (upon thermal averag- 
ing). As such, it suffices to consider a first-order per- 
turbation theory in imaginary time for the total density 
in the polarqn picture. With some details elaborated in 
Appendices [13, we finally obtain p]? « Pg^^g) -I- Ps^(i)i 
with 
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tanh(/377/2); (15) 

dTSn{T)K„{T), (16) 



'e^ + i\\. Here Snij) = {an{T)(J-) and 
{Vn{T)D'^_) Hb are the correlation functions of 
the system and of the bath, respectively. They are given 
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V ^ V 

X (ecosh[i(/3 - 2t)77] + ^ sinh[i(/3 - 2t)tj]^ , 

(17) 

-^sech(/377/2) 

X f cosh[i (/3 - 2r)7?] + - sinh[i(/3 ~ 2t)tj] 
\ 2 n 2 



(18) 

The bath correlation functions are Kx{t) = 2Cx{t)/A 
and Ky{T) = 2zCy(r)/A. Note that the first-order cor- 
rection here is again linked with the above-defined bath 
correlation function C„(t). So by construction, our per- 
turbation theory for off-diagonal elements of RDM works 
even better for stronger SBC coupling (thus analogous to 
the previous treatment for and p'^). 



III. RESULTS 

With all the RDM elements analytically obtained 
above, we now examine the validity of our theoretical 
results and reveal interesting physics in non-canonical 
equilibrium statistics. Instead of examining all the RDM 
elements (one exception later), we use a single quantity 
to characterize non-canonical statistics, i.e., the small- 
est possible angle (9) to be rotated (in radians) on the 
Bloch sphere to reach the eigenstates of Hs from the 
diagonal representation of the RDM. The theoretical re- 
sults (solid line) are plotted in Fig. [T] as a function of 
the SBC strength, 7, for a fixed temperature. For small 
values of 7, 9 is small, so the RDM's diagonal repre- 
sentation is close to, or significantly overlapping with, 
that of Hs- This is expected because for weak SBC 
strength, the equilibrium statistics should be canonical. 
As 7 increases, 9 increases, indicating that the RDM di- 
agonal representation continuously and monotonously ro- 
tates away from the eigenstates of Hs [iJl- To further 
elucidate the continuous change in 9, an analogous an- 
gle, namely, the angle the RDM diagonal representation 
should be rotated to reach the eigenstates of Hs b is also 
plotted in Fig. [T] (dashed line). Interestingly, for large 
values of 7, the RDM diagonal representation is seen 
to approach that of Hsb- Indeed, because C„(r) —J- 
for large 7, all the pcrturbative corrections to pg^, p|?, 
and p]? approach zero, and hence the RDM approaches 
exp(— I/Jctj), whose diagonal representation should be 
parallel to that of Hsb (as both are a function of (Tz). 
This is the case at arbitrary temperature. Our theoret- 
ical results are also in quantitative agreement with the 
solid dots shown in Fig. [1] obtained numerically from 
Monte Carlo simulations based on imaginary time path 
integral (a powerful method if the bath temperature is 
not too low [1^). That is, for a varying SBC strength. 
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either weak or strong, our theory and numericaUy exact 
results agree. This confirms that our analytical treat- 
ment for RDM off-diagonal elements performs equally 
well in the regime valid for treating the RDM dia gon al 
elements (hence almost the entire parameter space |1C| ) . 



1.2 




0.1 0.2 
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FIG. 1: Coupling-strength dependence of the angle to be ro- 
tated on the Bloch sphere to reach eigenstates of Hs (solid 
line) or Hsb (dashed line) from eigenstates of equiUbrium 
RDM, for /3 = 1, e = 0.5 and tj^ = 5 (in unit of A). Sohd 
dots are numerically exact Monte Carlo simulations results 
(for details of this method, see Ref. [11). 

In a different context, i.e., decoherence dynamics 
[l6l [TtI . the RDM diagonal representation is regarded 
as a special representation, often called a preferred basis 
of decoherence. It is in this special representation that 
decoherence can be understood as the disappearance of 
the off-diagonal matrix elements of a time-evolving RDM. 
A recent study using a low-dimensional quantum chaos 
model as a quantum bath [3l shows that the preferred 
basis of decoherence shows exactly the same qualitative 
behavior as observed in Fig. [Ijfor the equilibrium RDM, 
i.e., after a short period of decoherence the preferred 
basis of a system coincides with the eigenstates of Hs 
for weak SBC and becomes the eigenstates of Hsb f^or 
strong SBC, with a continuous deformation in intermedi- 
ate regimes. This very feature shared by the equilibrium 
RDM considered here and the preferred basis of decoher- 
ence is somewhat expected: an equilibrium RDM is an 
asymptotic result of quantum dissipation. Due to this 
interesting connection, the particular diagonal represen- 
tations of RDM as a result of non-canonical statistics 
can be also understood as one remarkable (previously 
overlooked) outcome of Nature's super-selection in open 
quantum systems [l^, [l9l - |22j . 

We now turn to the temperature dependence of non- 
canonical statistics at a fixed 7, as depicted in Fig. [2] for 
7 = 0.1. We choose 7 = 0.1 as an example because it 
represents an intermediate SBC strength in Fig.[TJ As ob- 
served from Fig. [5J for temperature lower than ksT = 1 
(a value considered in Fig. [1]) , the RDM diagonal repre- 
sentation is further rotated from that of Hs (solid line) , 
but gets closer to that of Hsb (dashed line). There- 
fore non-canonical statistics becomes more pronounced 
when temperature decreases. On the other hand, when 
temperature increases, the plotted angles continuously 



change in opposite directions, showing that the RDM 
diagonal representation gradually moves away from the 
eigen-representation of Hsb but smoothly approaches 
that of Hs- For temperature values much higher than 
shown in Fig. [21 this trend persists. Numerically exact 
Monte-Carlo simulation results (solid dots) are also pre- 
sented in Fig. m thus supporting again our theory. 




FIG. 2: Temperature dependence of the angle to be rotated 
on the Bloch sphere to reach eigenstates of Hs (solid line) or 
Hsb (dashed line) from eigenstates of equilibrium RDM, for 
7 = 0.1, e = 0.5 and lOc = ^ (in unit of A). Sohd dots are 
numerically exact Monte Carlo simulations results. 

Theoretical results outlined above may be further ex- 
ploited to understand the asymptotic high-temperature 
behavior of RDM. Keeping only terms of /3° and /?^, 
the use of Eq. ^ gives p^^ = 1/2 + and pf = 
1/2-6/3/4; whereas the use of Eqs. ([T5)) and (fT6)) yields 
p^s = —A/3/4 [l3|. But these asymptotic density matrix 
elements are exactly those of a canonical distribution for 
Hs = |crz + yCTx at high temperature. Non-canonical 
statistics is seen to have totally vanished at high temper- 
ature. This clearly depicts an interplay between two com- 
peting factors: SBC strength and temperature. A larger 
SBC strength generates a non-canonical RDM, but the 
thermal averaging tends to dilute non-canonical statistics 
and wipes it out completely at high temperature. Inter- 
estingly, this competition can be also appreciated from 
the bath statistics which is also non-canonical. In partic- 
ular, in terms of the boson occupation number on mode fc, 
the ratio of the leading-order correction to the canonical 
result is proportional to gf/ (ukkBT) p^ . which becomes 
negligible at high temperature. 

To further examine non-canonical statistics character- 
ized by a single angle measure (6), we present in Fig. O 
d9/dT, i.e., the sensitivity of the RDM diagonal repre- 
sentation to a small temperature variation, as a function 
of temperature. The sensitivity is low for very low tem- 
perature, but it rapidly increases, reaching a maximum 
at characteristic temperature scales that are comparable 
to other system parameters (such as 7). As temperature 
increases further, the sensitivity drops to zero asymptot- 
ically. The sensitivity profile 0jS 0j function of temper- 
ature qualitatively changes for a varying SBC strength. 
For weak SBC strength (e.g., 7 = 0.05), it exhibits a 
sharp peak. Hence the rotation of the RDM diagonal 
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representation mainly occurs within a narrow tempera- 
ture window. For strong SBC strength (e.g., 7 = 0.2), 
the sensitivity profile displays a rather flat structure, sug- 
gesting that it is harder for thermal effects to compete 
with SBC. Thus, if and only if the SBC strength is rather 
weak, then the temperature that gives cPO/dT'^ = (i.e., 
largest sensitivity of the RDM diagonal representation to 
temperature) becomes an interesting temperature scale. 




FIG. 3: Sensitivity of the RDM diagonal representation to 
temperature variation, as described by dO/dT vs T. is the 
angle to be rotated to reach eigenstates of Hs from eigenstates 
of the equilibrium RDM. System parameters are given by e = 
0.5 and lOc = 5 (in unit of A). 



IV. CONCLUSIONS 

For an open quantum system not weakly coupled with 
a bath, its equilibrium state is far from a canonical state 
at low temperature. Exact analytical solutions are typ- 
ically not available (one known exception is the model 
of a harmonic oscillator linearly interacting with a boson 
bath). A systematic approach to such open quantum sys- 
tems is hence highly desirable in efforts to better under- 
stand their qualitative and quantitative features of equi- 
librium statistics as temperature and/or SBC strength 
varies. To our knowledge, the theoretical treatment ad- 
vocated in this work, as supported by numerical results, 
represents the first attempt along this direction that can 
almost cover the whole range of SBC strength and the 
whole range of temperature of a spin-boson model. Be- 
cause non-canonical statistics is closely related to strong 
system-bath correlation, we anticipate our theory to be 
also useful in understanding system-bath entanglement. 

Our theoretical findings based on a spin-boson model 
are much relevant to experiments based on quantum dots. 
Acoustic phonon modes have been identified as the prin- 
cipal source of decoherence in InGaAs/GaAs quantum 
dots [131 , and temperature is widely tunable in such 
a semiconductor implementation. Certainly, in a real 
system the bath spectral density may not be the super- 
ohmic one assumed here. To address this concern we have 
carried out numerically exact calculations for an Ohmic 
bath at a nonzero temperature, obtaining results that are 



qualitatively the same as presented in this work [13| (even 
though our analytical treatment based on a full polaron 
transformation cannot be applied to this case). Our ap- 
proach can also be generalized to study any dissipative 
large spins as well as to ensembles of spins coupled to a 
common bath 1251. 
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Appendix A: Derivations of the Off-diagonal RDM 
Element 



The off-diagonal element of the equilibrium reduced 
density matrix (RDM) can be formally written as = 

tr[cr_pf] = tT[a^D'ipt] where pt = tJ^tStt- Expanding 



tr[e-/3«t] ■ 

Pt up to first order in Hsb, we have 



where 



Pt ~ Pt,{Q) + Pt,(i), 



-pHo 



pt.il) 



(Al) 



(A2) 



-l3Ho rP 



tr[e-/5«o] 7o 



dTHsBir). 



(A3) 



Inserting the above expressions into p^^? = tr[a^D^pt], 
wc obtain 



2^ 



■tanh(/3?7/2); 



(A4) 



Pf,(l) = tr[cr_i:'^Pf,(i)] 



n—x.y ^ 

= - E dr{an{r)a.)^jVr,{T)D'_)Hs, 
= - I dTSnir)K^{r), (A5) 

n—x,y ^ 

where the explicit expressions for the system and bath 
correlation functions, 5'„(r) and Kn{T), are already given 
in the main text. 
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Appendix B: High- Temperature behavior of RDM 

Here we give some details to see how non-canonical 
statistics of RDM totally vanishes at high temperature. 
For the off-diagonal element, the zeroth order term van- 
ishes at high temperature since R decays exponentially 
with T. Thus, only the first-order correction term, p^f^iy 
contributes to the off-diagonal element at high temper- 
ature. Furthermore, at high temperature the system 
correlation functions can be approximated as Sx{t) sa 

lj,(/3-2T)e/2 g^^^ ^ „ilc(/3-2r)e/2 (^^CrC WC havC 

used r] ~ e). Note also that though R vanishes at high 
temperature, the term i?^e'^^'^' contained in the bath cor- 
relation functions remains finite (keeping in mind that 
< T < That is. 



exp 



duj J{uj) 



cosh[i/3w] - cosh[i(/3 - 2t)lu] 



exp 



sinh(/?w/2) 

r2 - I3t1 



(Bl) 



This remarkable agreement nicely demonstrates that 
canonical statistics is recovered at high temperature. 



Appendix C: Ohmic Bath 

Here we study the rotation angle between the eigen- 
states of the RDM and the eigenstates of iJg or Hsb 
using an Ohmic bath, J {to) = 7a;c~"/"'=. Unfortunately, 
a full polaron method as we used in the main text is 
not applicable for an ohmic bath as it suffers from an 
unphysical divergence issue 0, [1^ . The integral in the 
renormalization constant R is divergent for all coupling 
strength, and the tunneling element is always normalized 
to zero. Therefore here we only present the numerical re- 
sults from the imaginary time path integral simulations 
(for not too low temperature). The coupling and tem- 
perature dependence of the rotation angle are plotted in 
Fig. |4] and Fig. [5l It can be seen that the features of the 
figures are qualitatively similar to those obtained using a 
super-ohmic spectral density in the main text. Therefore, 
the general observations made in our main text should 
not be sensitive to the spectral density of the bath. 



where k = 4 = ^uj^ and we have used the 

expansions sinh(x) ~ x and cosh(a;) « 1 + ^a;^ to arrive 
at the second expression. The bath correlation functions 
at high temperature can then be written as Kx{t) « 
f exp[-K^^^] &ndKy{T) « iAexp[-Kl^]. We then 

have p'i « _A J^/3g{/3-2r)e/2-K(r^-/^r)//J^^^ gj^^^g ^ ig 

also small (0 < t < /3), the integrand can be further 
expanded using k, 1 + x and we finally obtain pg' up 
to the second order in fi: 
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Ps 



A 



(/3 - nP^/Q). 



(B2) 



Calculating the diagonal elements at high temperature 
is more straightforward: the double integral in matrix A 
[see Eq. (12)] indicates that it is at least proportional to 

and can be discarded if we are only interested in terms 
up to the first-order of /3. The diagonal elements, p]^ and 

can then be written as ^(1 ± It: s[cr zPs'^\) ■ The are 
explicitly given by 



Ps 



i[l--tanh(/?r;/2)]«i(l-e/?/2); (B3) 
2 T] 2 

l[l + -tanh(/?r;/2)] « hl + e(3/2), (B4) 
2 T] 2 



where we have used 77 = -^/e^ + Ajj « e and tanh(a:;) w x. 

Gathering all the results above, our analytic theory 
predicts that at high temperature. 
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FIG. 4: The angle to be rotated on the Bloch sphere to reach 
the eigenstates of //s (crosses) or Hsb (solid dots) from the 
eigenstates of the equihbrium RDM as a function of coupling 
strength, for 13 — 1, e = 0.5 and ojc = 5 (in unit of A). 
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The above expression is exactly the same as the high- 
temperature canonical state of the system, 1^. . . 



FIG. 5: The angle to be rotated in the Bloch sphere to reach 
the eigenstates of (crosses) or Hsb (solid dots) from the 
eigenstates of the equilibrium RDM as a function of temper- 
ature, for 7 = 1.5, e = 0.5 and lOc = 5 (in unit of A). 



7 



Appendix D: Eigenvalue at Zero Temperature 



Appendix E: Bath Statistics 



In the eigcnbasis, the equihbriuni RDM can be written 
as PS = {'^^ X2) ^ where are the eigenvalues of RDM 
and Ai + A2 = 1. The eigenvalues denote the population 
of each of the eigenstate. The eigenvalues also serve as an 
indicator of the purity of the system. If both eigenvalues 
are non-zero, the system is in a mixed state. Below, we 
will use the larger eigenvalue, A2, to investigate the purity 
of the system: the system is in a pure state if A2 = 1 and 
vice versa. 

Due to the finite system-bath coupling, the equilibrium 
RDM might not be a pure state at T = even though the 
system plus the bath is in their entangled ground state. 
To examine how the purity of RDM at T = depends 
on the coupling strength, we plot A2 as a function of 7 
in Fig. [HI It can be observed that the eigenvalue exhibits 
an interesting non-monotonic behavior as a function of 
the system-environment coupling strength. At 7 = 0, 
the system is in the ground state of Hs with unit pu- 
rity. At finite coupling, both eigenstates are populated 
and RDM is a statistical mixture due to the system-bath 
entanglement. Interestingly, RDM is reduced to a pure 
state at very large 7, indicating that the system-bath en- 
tanglement vanishes at ultra-strong system-environment 
coupling. However this pure state is no longer the eigen- 
state of the system Hamiltonian, but that of in the 
interaction Hamiltonian. 



Here we examine the equilibrium statistics of the bath 
by examining the average boson number of each mode, 
which is denoted by (n^). In the polaron frame, the boson 
number operator is given by 



rik 



F 

e rifeC 



9k 



9l 



<^M + h) + ^, (El) 

where F = az X]fe ^(^1 ~ ^1^) ^^'^ ^^^^ line. An ap- 
proximate expression of (rik) can be obtained by 



ink) 



tr e 



-PHt] 



trfe^rifeC ^ e^e ^^*e 



trfn/cC 



tr[e-/3ff*] 



tr[e 
trfrifcC 



tr[e-/3«o] 



(E2) 




FIG. 6: The larger eigenvalue, A2 of RDM, plotted against 
the SBC strength 7, for T — 0, e = 0.5 and uJc ~ 5 (in unit of 
A). 



Inserting Eqn. (|E1[) into the above expression, we have 



(E3) 



(nk) = {nk)o + gl/ujt 



where (n/j)o = (e'^s^ — 1)^^ is the average boson num- 
ber without system-bath coupling. In the high tem- 
perature limit, it satisfies the equipartition theorem, 
('T-fc)o ~ ^isZ Therefore, the fractional correction, given 



'{nk)o 



perature. 



9k 



ujkki 



7p, becomes negligible at high tem- 



[1] T. Brixner, J. Stenger, H. M. Vaswani, M. Cho, R. E. 

Blankenship, and G. R. Fleming, Nature 434, 625 (2005). 
[2] J. Q. You and F. Nori, Nature 474, 589 (2011). 
[3] See, for example, K. Hennessy, A. Badolato, M. Winger, 

D. Gerace, M. Atatiire, S. Guide, S. Fait, E. L. Hu, 

and A. Imamoglu, Nature (London) 445, 896 (2007), G. 

Giinter, A. A. Anappara, J. Hees, A. Sell, G. Biasiol, L. 

Sorba, S. De Liberato, C. Ciuti, A. Tredicucci, A. Leit- 

enstorfer, and R. Huber, Nature 458, 178 (2009) and A. 

Auer and G. Burkard (2011), arXiv:1112.3211. 
[4] One remarkable example is specific heat anomalies, see, 

for example, P. Hanggi, G. L. Ingold, and P. Talkner, 



New J. Phys. 10, 115008 (2008) and G. L. Ingold, 
P. Hanggi, and P. Talkner, Phys. Rev. E 79, 061105 
(2009). 

[5] Technically, this is because in the classical limit of a wide 
class of microscopic open-system models, the Hamilto- 
nian of mean force coincides with the bare system Hamil- 
tonian (which is generically not the case in open quantum 
systems). See, for example, M. Campisi, P. Talkner, and 
P. Hanggi, Phys. Rev. Lett. 102, 210401 (2009) and M. F. 
Gelin and M. Thoss, Phys. Rev. E 79, 051121 (2009). 

[6] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 2008). 



8 



[7] H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 

168, 115 (1988). 
[8] D. Braak, Phys. Rev. Lett. 107, 100401 (2011). 

[9] We use {■)hb to denote an average over ^j^jg-A^^^j • 
Similarly, later we use (■)^^ denoting an average over 

trs[e-''-Ss] ■ 

[10] C. K. Lee, J. Moix and J. Cao, J. Cliem. Phys. 136, 
204120 (2012). 

[11] This expression without the polaron transformation was 
previously used by B. B. Laird, J. Budimir, and J. L. 
Skinner, J. Chem. Phys. 94, 4391 (1991) and E. Geva, 
E. Rosenman, and D. Tannor, J. Chem. Phys. 113, 1380 
(2000). 

[12] Using a super-ohmic spectral density, we have (/)(r) = 

[13] See Appendix for necessary details of our theoretical 
derivations, high-temperature behavior of RDM, addi- 
tional numerical results for a different bath spectrum 
density, dicussions on the eigenvalues of RDM as a func- 
tion of system parameters, as well as the bath statistics. 

[14] For a spin-boson model with sub-ohmic or ohmic spectral 
density, one may expect an abrupt phase transition at 
very low temperature. 

[15] J. M. Moix, Y. Zhao, and J. Cao, Phys. Rev. B 85, 



115412 (2012). 
[16] W. H. Zurek, Phys. Rev. D 24, 1516 (1981). 
[17] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 
[18] W.-G. Wang, L. He, and J. B. Gong, Phys. Rev. Lett. 

108, 070403 (2012). 
[19] J. P. Paz and W. H. Zurek, Phys. Rev. Lett. 82, 5181 

(1999). 

[20] D. Braun, F. Haake, and W. T. Strunz, Phys. Rev. Lett. 

86, 2913 (2001). 
[21] W.-G. Wang, J. B. Gong, G. Casati, and B. Li, Phys. 

Rev. A 77, 012108 (2008). 
[22] C. Gogolin, Phys. Rev. E 81, 051127 (2010). 
[23] A. J. Ramsay, A. V. Copal, E. M. Ganger, A. Nazir, 

B. W. Lovett, A. M. Fox, and M. S. Skolnick, Phys. Rev. 

Lett. 104, 017402 (2010). 
[24] A. J. Ramsay, T. M. Godden, S. J. Boyle, E. M. Ganger, 

A. Nazir, B. W. Lovett, A. M. Fox, and M. S. Skolnick, 

Phys. Rev. Lett. 105, 177402 (2010). 
[25] See, for example, T. Vorrath and T. Brandes, Phys. Rev. 

Lett. 95, 070402 (2005) and L. D. Contreras-PuUdo and 

R. Aguado, Phys. Rev. B 77, 155420 (2008). 
[26] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 

Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 

1 (1987). 



